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ABSTRACT 

The recently discovered Galactic X-ray transient XTE J 1752— 223 entered its first 
known outburst in 2010, emitting from the X-ray to the radio regimes. Its general X- 
ray properties were consistent with those of a black hole candidate in various spectral 
states, when ejection of jet components is expected. To verify this, we carried out very 
long baseline interferometry (VLBI) observations. The measurements were carried out 
with the European VLBI Network (EVN) and the Very Long Baseline Array (VLBA) 
at four epochs in 2010 February. The images at the hrst three epochs show a moving 
jet component that is significantly decelerated by the last epoch, when a new jet 
component appears that is likely to be associated with the receding jet side. The 
overall picture is consistent with an initially mildly relativistic jet, interacting with 
the interstellar medium or with swept-up material along the jet. The brightening of the 
receding ejecta at the final epoch can be well explained by initial Doppler deboosting 
of the emission in the decelerating jet. 

Key words: stars: individual: XTE J 1752—223 - stars: variable: others - ISM: jets 
and outflows - radio continuum: stars - X-rays: binaries. 
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1 INTRODUCTION 

It is clear from the literature of the past three decades 
that, for almost every black hole X-ray transient (XRT) 
observed at radio wavelengths , a radio counterpart has 
been discovered (jFenderl 120060 . In a small number of 
sources, the ejecta have been resolved and monitored 
as they travel away from the central source. Thus, it 
has been possible on rare occasions to measure proper 
motions, sometimes at appa rently super-luminal veloc- 
ities: e.g. GRS 1915+105 (|Mirabel feFtodriguezl Il994j ; 



Fender et all 119990 GRO J1655-40 (Tingav et all Il995l: 
Hiellming fc Rupenl Il995l ). GX 339^4 l|Corbel et all 



2010h . and even detect the only known instances of a 
Galactic parsec-scale jet being de c elerate d in the ISM: 
XTE 
2003h 



J1550-564 
and XTE 



(|Corbel et alj 120021: iKaaret et all 



Miller- Jones et al.l 



J17 48- 
2008). 
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Jet 



l|Hiellming et all 1 1998t 
deceleration has been 



also investigated in GRS 1915+105 although no conclusive 
evidence was found J Miller- Jones et aT1l2007| ). 
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The X-ray transient XTE J1752-223 was discovered by 
the Rossi X-ray Timing Ex plorer (RXTE) on 2009 October 
23 l|Markwardt et all 120091 ') at the start of its first known 
X-ray outburst. It showed a long and gradual rise at X-ray 
energies, whilst remaining sp ectrally hard. The X-ray source 
later e volved, became softer {Homan 20101: iBrocksopp et all 
l2010al ) and entered a spe ctral state common l y asso ciated 
with jet ejection events l|Fender et all |2004 120091 ). The 
outburst has been well mo nitored by the Monito r of All- 
sky X-ray Image (MAXI , iNakahira et~al1 l201Ch. RXTE 
jShaposhnikov et al.ll20~Lcl ), and Swift jCurran et afll2010l ) 
at X-ray energies. All these X-ray observations show that 
XTE J1752-223 is likely to be a black hole transient. 

Following the activation of XTE J1752-223, we initi- 
ated the ATCA (Australia Telescope Compact Array) radio 
observations and discovered the ra dio counterpart with a 
flux density of ~ 2 mjy at 5.5 GHz (|Brocksopp et afl [2009). 
The later ATCA monitoring observations showed that the 
radio source entered a series of flares, peaking at 5 - 20 
mjy (Brocksopp et al. in prep.), after the transition from 
X-ray hard to soft state. To detect the potential ejecta and 
study their evolution, we carried out an EVN (European 
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VLBI Network) experiment and three follow-up VLBA ex- 
periments at 5 GHz in 2010 February. In this paper, we 
present the results of these VLBI (Very Long Baseline In- 
terferometry) observations. 



2 RAPID EVN AND VLBA OBSERVATIONS 

The performed VLBI experiments are summarised in Ta- 
ble Q] The low declination and potentially weak flux of 
XTE J1752-223 were possible problems for VLBI obser- 
vations. Therefore, to quickly resolve these concerns, we 
initiated an e- VLBI experime nt with the western Euro- 
pean telescopes (|Szomorull2008l ) . The EVN experiment used 
1024 Mbps data rate (16 channels, 16 MHz per chan- 
nel, 2 bit sampling, dual polarisation). The real-time cor- 
relation was done with the Earth orientation parameters 
(EOP) predicted from the EOP model of one day earlier. 
We applied 2-second integration time, and 16 frequency 
points per channel. The participating stations were Medic- 
ina, Yebes, Torun, Onsala, and Westerbork. In the EVN 
experiment, we used the following J2000 coordinate: RA = 
17 h 52 m 15!095, DEC = -22°20'32'.'782 (positional uncer- 
tainty: a = 0"3), which was deter mined by the ATCA ob- 
servations (|Brocksopp et al J 120091 ) . To obtain fringe-fitting 
solutions and a reference point, we scheduled a nearby 
(0?8) phase-referencing source: PMN J1755-2232 (RA = 
17 h 55 m 26=285, DEC = -22°32'10'.'593, J2000, position er- 
ror a ~ 15 mas). As these sources have low elevation (< 30°) 
in Europe, we used a short cycle time: 160 seconds on the 
target and 80 seconds on the reference source. We also ob- 
served a strong and compact source (NRAO 530) as the 
fringe finder and bandpass calibrator. 

We successfully detected a radio sourc e, consistent with 
an eje ction from the black hole candidate (|Brocksopp et all 
l2010bT ). during the EVN experiment and then performed 
three follow-up VLBA experiments. We used the same 
phase-referencing source, cycle time, and observing fre- 
quency. The recording data rate was 512 Mbps (16 channels, 
8 MHz per channel, 2 bit sampling, dual polarisation) . There 
were eight VLBA telescopes available at the first epoch, nine 
at other two epochs. The data were correlated with the same 
parameters as the EVN experiment. 

We also performed an EVN experiment in March and 
two VLBA experiments in April to study other ejection 
events and attempt to detect the core. These additional 
results will be presented by Brocksopp et al. in a general 
paper. 



3 VLBI DATA CALIBRATION 

We used the NRAO software AIPS (Astronomical Image 
Processing System) to perform the initial calibrations. The 
a-priori amplitude calibration was done with measured sys- 
tem temperatures and antenna gain curves. We corrected 
the EOP model for the VLBA data before any phase cali- 
brations. We follow the same procedure for both the EVN 
and VLBA data reduction. (1) We corrected the parallactic 
angle. (2) We ran the global fringe-fitting for NRAO 530 
with half-minute solution interval and then solved for the 
instrumental bandpass. (3) With the bandpass solutions, 



Table 1. The summary of the image parameters of Fig. [T] 
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number of the used telescopes; (5) total observing time; (6) peak flux density; 
(7) off-source noise level; (8 - 10) sizes of the major and minor axes of restoring 
beam and its position angle. 

we re-ran the fringe-fitting to solve for the instrumental 
phase and delay using two-minute data of NRAO 530. (4) 
We ran the fringe-fitting to solve for the phase, the fringe 
rate, and the delay for the calibrators with a solution inter- 
val of the scan length (~ 1 minute). On the long baselines 
to Mauna Kea (Mk) and St. Croix (Sc), there were no so- 
lutions found for PMN J1755— 2232 as it has much lower 
correlation amplitude (< 10 mjy) most likely due to scat- 
ter broadening. In view of this problem, we removed both Sc 
and Mk. The phase wrapped slowly (fringe rate < 5 milliHz). 
The solutions of PMN J1755-2232 were then transferred to 
XTE J1752-223 by linear interpolation. (5) The data were 
averaged in each IF and then split into single-source files 
after all the corrections were applied. 

The reference source PMN J1755— 2232 was imaged 
by circu lar Gaussian model fitting and self-calibration in 
DlFMAP( |Shepherd et al.l Tl994L The source was well repre- 
sented by a circular Gaussian model with a size of 4.2 mas 
and a total flux density of ~ 200 mjy at 5 GHz. Finally, we 
self-calibrated the u-v data with the model and applied the 
amplitude and phase solutions to both sources in AIPS. 



4 VLBI DETECTION OF XTE J1752 223 

The imaging results for the X-ray transient XTE J1752-223 
are shown in the left panel of Fig. [T] The top small pan- 
els from left to right are the EVN image of 2010 Febru- 
ary 11 and the VLBA images of 2010 February 18 and 26. 
The background large image is the VLBA image of 2010 
February 26. The cross in each image marks the position 
of component A at the first epoch (RA = 17 h 52 m 15 ! !06370, 
DEC = -22°20'31'.'9838, J2000). At the bottom of the large 
background image, there is another jet component marked 
as component B at an angular separation 488 mas from com- 
ponent A and at a position angle consistent with the moving 
direction of component A. The related map parameters are 
listed in Tab. [1] 

Component A, surrounded with the beam pattern, is 
clearly seen with a peak brightness of 2.32 mjy beam -1 in 
the dirty map at the first epoch, when natural weighting 
is used. After removing component A with a circular Gaus- 
sian model, we notice that there may be at least one more jet 
component. One candidate is located at angular separation 
18.7 mas, position angle — 84? 0; the other at angular sepa- 
ration 70.6 mas and position angle —36? 3. Both candidates 
have a peak brightness ~ 0.91 mjy beam -1 (~ 4er rms ) using 
natural weighting. In Fig. []J,, the two candidates show the 
second positive contours. If either candidate is removed by 
circular Gaussian model fitting, the other also becomes faint. 
If we reduce the contribution of the long baselines, both be- 
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Figure 1. The decelerating jet of the X-ray transient XTE J1752— 233. All the VLBI images are centred at the location of component A 
on 2010 February 11, indicated by a cross. Component B was not detected until 2010 February 26. The contours start from 3<r off-source 
noise level and increase by a factor of -1.4, -1, 1, 1.4, 2, 2.8. The related map parameters are listed in Tab. IT1 The right panel plots the 
fitting results (Tab. [3j of the proper motion data (Tab. [2j °f component A using the models with (solid curve) and without deceleration 
rate (dotted line). The reference time MJD 55238.4 corresponds to the first VLBI observations. 



come brighter and show a small peak (~ 5%) brightness 
difference. Due to the limited sensitivity and u — v coverage 
during the 1.2-hour observations, neither components can be 
unambiguously identified as a true jet component. However, 
there is evidence for the extended emission for the source 
as the total restored flux density is much lower than that 
(~ 16 mjy at 5.5 GHz) measured by the ATCA (Brocksopp 
et al. in prep.). 

In the follow-up VLBA observations, the higher reso- 
lution and sensitivity are achieved by more telescopes and 
longer observing time. To image the extended source, we 
used natural weighting again. Because of the resolved struc- 
ture and the decaying peak flux density, the source is only 
seen clearly in the dirty map with the synthesised beam 
(16.2 x 3.9 mas at position angle — 15?6) at the first VLBA 
epoch. However, the large-scale beam pattern around the 
faint source could also be easily identified at the later two 
epochs. If we taper the long baselines, use the short baselines 
only, or increase the image pixel size, the source becomes sig- 
nificantly brighter in the dirty map at the later two epochs. 
None of the suspected ejecta candidates in the EVN image 
are further seen after 7 days in the later VLBA images. Be- 
cause the diffuse emission can not be well restored by clean 
components, Gaussian models were used in making all the 
VLBI images of Fig. Q] Due to the limited SNR (6 - 12), 
circular rather than elliptical Gaussian model fitting was 
adopted to reduce the number of free parameters. 

Table [2] lists the best-fit parameters of the circular 
Gaussian model. To show the motion of component A, we 
take the position of component A measured at the first epoch 
as the reference origin. The random position error was es- 
timated by ^ 2Snr" ~ ' wnere ^ ma l an d bmin are the size of 
the major and minor axes of the used restoring beam, SNR 
is the signal to noise ratio (^2*) listed in Column (7) of 
Table [2] Note that the rather large systematic position error 



Table 2. The circular Gaussian model fitting results of the de- 
tected jet components in the X-ray transient XTE J1752— 223. 
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from the reference source will not affect our proper motion 
measurements. The fitted size has the same random error as 
the position for each component. Since the measured sizes 
(^ 8 mas) are much larger than that (4.2 mas) of the ref- 
erence source, they should be very close to the true size of 
the ejecta. At the second epoch, we notice that component 
A shows an elongated structure in the East- West direction 
and the eastern side is brighter than the western side. This 
brightness distribution caused a slightly different position 
angle of the component, compared to what is measured at 
later epochs. The VLBI flux density errors are usually ~ 5%. 



5 GRADUAL JET DECELERATION 

The angular separation of component A versus time is shown 
in the right panel of Fig. [1] We take the position and the 
time of component A measured at the first epoch as the 
reference origin. We fit these data points to the following 
proper motion model: 

r = ro + not — 0.5/it 2 (1) 

where r is the angular separation; t is the observing time; 
ro and fio are the angular separation and the proper mo- 
tion at t = 0, fi is the apparent deceleration rate. The dot- 
ted straight line and the solid curve represent the uniform 
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Table 3. Best-fit parameters using the proper motion models 
with and without the deceleration rate. 
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proper motion model (fi = 0) and the proper motion model 
with the deceleration rate (fj, =^= 0) respectively. The best-fit 
parameters are listed in Table [3] The model of ji = gives 
an average proper motion of fl = liq = 6.90±0.05 mas day -1 
with the reduced \ 2 — 118.4. The degree of freedom (DoF) 
was listed in the last column. The model of fx ^ gives 
= 9.15 ± 0.15 mas day -1 at MJD 55238.4 and a deceler- 
ation rate of fi = 0.34 ± 0.02 mas day -2 with the reduced 
X 2 = 3.9. With the deceleration rate, component A has a 
proper motion of 4.0 mas day -1 at the last epoch. It is clear 
that the deceleration rate should be included in the proper 
motion model. 

Jet deceleration was a lso found in XTE J1550— 564 us- 
ing C handra observations l|Corbel et al.|[2002l ; iKaaret et al.l 
l2003h and XTE Jl 748-288 with VLA observations 
l|Hiellming et all 1 19981 ; iMiller- Jones et all 120081 '). Compared 
with them, the observed deceleration in XTE J1752— 223 is 
free from the blendi ng of multiple .jet compon ents caused by 
the low resolution (|Hiellming fc Rupenlll995r ). If there is a 
sequence of ejecta which decreased sequentially more rapidly 
in flux density with increasing distance from the core, the 
cluster of components may show a decreasing proper mo- 
tion. In our case, these VLBI observations have a resolu- 
tion of < 10 mas and can well identify single ballistic ejecta 
with a time resolution of less than one day, assuming an 
initial proper motion 10 mas day . XTE J1752— 223 is the 
second known case of gradual jet deceleration, although on 
a much smaller scale (~ 100 mas) than the first case of 
XTE J1550— 564. As for the previous two sources, the jet 
deceleration in XTE J1752— 223 is most likely due to inter- 
action with the external dense interstellar medium (ISM) or 
the residual slowly-moving ejecta from the previous ejection 
along the jet path. 



6 COMPONENT B: THE RECEDING EJECTA? 

We interpret component A as an approaching jet component 
since it is the only component detected at the four epochs 
and our VLBI observations were performed just after the as- 
sociated radio flare reached its peak flux density. The ATCA 
observations (Brocksopp et al., in prep.) show that it had a 
decaying flux density and a pretty stable and steep spectral 
index: a — —1 (S v oc v a ) between 5.5 and 9 GHz, i.e. there 
was no indication of another ejection event during our VLBI 
observations. Thus, the possibility that the components A 
and B detected at the last three epochs are associated with a 
different ejection event can be ruled out. Besides the proper 
motion, component A shows a hint of expansion. The evo- 
lution of its size is displayed in Fig. [2] The first three data 
points give an expansion speed of 0.9 ± 0.1 mas day - with 
reduced \ 2 = LI- The expansion speed is much slower than 
its proper motion, indicating that its expansion was signifi- 
cantly confined too. Because component A shows an increas- 



ing size and a decaying peak brightness, its size estimation 
at a later stage is limited by the image sensitivity and the 
lack of short baselines. For this reason, the last data point 
was omitted in the linear fitting. According to the evolution 
of component A, a jet component is expected to have more 
compact structure at the earlier stage. Compared with the 
size (7.9 mas) of component A measured at the first epoch, 
component B shows a much larger size (11.9 mas). This in- 
dicates that component B is most likely an evolved compo- 
nent, which was ejected on the receding jet side at the same 
time as component A. 

According to the expansion speed, component A was 
ejected 8.7 days earlier, i.e. at MJD 55229.7 (2010 Feb. 2), 
which is represented by the x-intercept in Fig. [2] As com- 
ponent A may have significantly larger expansion speed and 
strong Doppler boosting effect if it is unhindered at this ear- 
lier stage, the inferred birth date may be the earliest possible 
birth date. Although such extrapolation is not guaranteed, 
the inferred date is at the beginning of the initial rising stage 
of the associated radio flare in the ATCA radio light curve 
(Brocksopp et al. in prep.). The average separation speed of 
the pair of components is papp + Mrcc ^ 20.4 mas day -1 
if they were indeed ejected on the inferred birth date. 
Since /I ap p ^ A*rcc, there is ft app > 10.2 mas day - , 
which is significantly higher than the average proper motion 
(6.9 mas day -1 ) measured during our observations. Thus, 
component A had already been significantly decelerated be- 
fore our VLBI observations. 

If the jet expansion is linear and symmetric on both 
sides, the ra tio of the approaching a nd receding component 
sizes is (e.g. IMiller- Jones et al.ll2004l ): 

-Rapp(tapp) _ tapp _ 1 + /3(0,^app) COS 9 
-Rrec(irec) t rC c 1 — /3 (0, t vcc ) COS 9 

where t apP and t TCC are the intrinsic times at which light 
leaves the approaching and receding jet components respec- 
tively and arrive at the telescope at the same observing 
time, P is the average jet speed in units of light speed c 
and 8 is the inclination angle of the jet axis. By a linear 
extrapolation, the approaching jet component has a size of 
21.3 mas at the fourth epoch. If there is no deceleration, 
/3(0,i ap p) = ,3(0, tree) = P, we can give /3cos# = 0.3c, which 
requires /3 ^ 0.3c and 9 < 73°. Note that the size of the re- 
ceding jet component detected for the first time is not likely 
to be affected by the over-resolution since it is much younger 
than the approaching jet component (t ICC ~ 0.5t apP )- Given 
the jet deceleration and t ICC ^ i apP at the same telescope 
time, then j3(0,t vcc ) > 0.3c and /3(0,i ap p) 5? 0.3c. Therefore, 
we can take 0.3c as the lower limit of the jet birth speed in 
the case of the jet deceleration. 

The ratio of the flux density measured at 7? a pp = 
Riac (iapp = troc, free from its intrinsic luminosity 
evolution effect) for a pa ir of discrete jet components 
l|Mirabel fc Rodriguez! [l99Sh : 

Sapp _ / 1 + P COS 8 \ 3 Q , . 

S Tec ~ v l-/3cos6> J ( ' 

The receding jet component had a size of R ICC = 11.9 mas 
at the last epoch. The corresponding observing time for the 
approaching jet component is at MJD 55242.9 (between the 
first two epochs). At t apP = t ICC , the radio core is inferred 
to be at the centre: angular separation ~ 174 mas and posi- 
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Figure 2. Evolution of the size of component A. The solid 
straight line shows the fitting result of the first three data points. 

tion angle ~ —130°. The radio core was not detected during 
any of the VLBI epochs but, since all four observations took 
place during the X-ray soft stat e, this is to be expec t ed, ac - 
cording to the unified model of iFender etafl (|2004 l2009h . 
The radio position con firms that its optica l counterpart is 
Swift-UVOT source A (|Curran et alJl201Ch . If we take the 
flux density at the first epoch as the upper limit of Sa. PP , 
then ^£ < 5 and /3 cos 69 < 0.2c < /3(0, t ICC ) cos 9, in agree- 
ment with the jet deceleration scenario on both sides. The 
observed flux density from the receding jet is deboosted by a 
factor: 8*£*, where 5 rcc = (1 - /3 2 )°' 5 (1 + /3 cos 0)" 1 . Because 
of the jet deceleration, the receding ejecta is less beamed 
away from our line of sight, and thus it looks relatively 
brighter. Note that the non-detection of the receding jet 
component at the first epoch may also be because it stayed 
at an earlier stage (t ICC ~ 0.5£ app if /3 cos 9 = 0.3) and its 
flux density was still much lower than the peak flux den- 
sity of the radio flare. The caveat in the above argument is 
that component B might have not followed the same lumi- 
nosity evolution model as component A. It is also possible 
that the brightening of the receding ejecta was due to sud- 
den jet-cloud interaction, as in the case of XTE J1748— 228 
(Hicllming ct al. I ll998i : iMiller- Jones et ai]|2008h . 

7 CONCLUSIONS 

In this paper, we present the results of the first VLBI 
observations of the new Galactic black hole candidate 
XTE J1752-223 during its first known outburst. With EVN 
and VLBA observations at four epochs in 2010 February, 
we imaged its radio counterpart at 5 GHz. We detect an 
ejected component at the first three epochs and find that 
its proper motion shows significant deviation from the uni- 
form proper motion model and requires a deceleration rate 
of 0.34 ± 0.02 mas day -2 . In the jet deceleration scenario, 
it has proper motion decelerating from 9.2 mas day -1 at 
the first epoch to 4.0 mas day -1 at the last epoch. It also 
shows slow but detectable variation of its transverse size 
indicating that its expansion is also significantly confined. 
This is the first time that a Galactic jet is found to be grad- 
ually decelerating on the hundred milliarcsecond scale. The 
discovery provides strong evidence for the existence of sig- 
nificant interaction around the jet at an early stage of its 
evolution. In addition to the approaching jet component, we 
detect another jet feature at the last epoch, which is most 
likely associated with the receding ejecta. We infer that the 
jet deceleration should start at a time much earlier than 
our VLBI observations using the birth date (around 2010 
February 2) from the ATCA radio light curve. Furthermore, 



we interpret the detection of the receding ejecta as a result 
of the reduced Doppler deboosting effect caused by the jet 
deceleration on the receding side and give a lower limit of 
0.3c for the jet birth sp eed assuming symmetric je t motion. 
It has been reported bv lShaposhnikov et alj (|2010h that the 
distance, estimated by the spectral-timing correlation scal- 
ing technique, is around 3.5 kpc. Thus, the proper motion 
observed at the first epoch would correspond to an apparent 
jet speed of ~ 0.2c, in agreement with our results (but note 
that the technique is very model dependent). 
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